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ABSTRACT 

Stellar parameters are not easily derived from M dwarf spectra, which are dominated by complex bands of diatomic and 
triatomic molecules and do not agree well with the individual line predictions of atmospheric models. M dwarf metallicities are 
therefore most commonly derived through less direct techniques. Several recent publications propose calibrations that provide 
the metallicity of an M dwarf from its Ks band absolute magnitude and its V — Ks color, but disagree at the ±0.1 dex level. 
We compared these calibrations using a sample of 23 M dwarfs, which we selected as wide (> 5 arcsec) companions of F-, G-, 
or K- dwarfs with metallicities measured on a homogeneous scale and which we require to have V band photometry measured 
to better than ~0.03 magnitude. We find that the Schlaufman & Laughlin (2010, A&A, 519, A105+) calibration has the lowest 
offsets and residuals against our sample, and used our improved statistics to marginally refine that calibration. With more 
strictly selected photometry than in previous studies, the dispersion around the calibration is well in excess of the [Fe/H] and 
photometric uncertainties. This suggests that the origin of the remaining dispersion is astrophysical rather than observational. 

Key words, stars: fundamental parameters - stars: binaries - general - stars: late type - stars: atmospheres - stars: planetary 
systems 



1. Introduction 

M dwarfs are the smallest and coldest stars of the main se- 
quence. Long lived and ubiquitous, M dwarfs are of inter- 
est in many astrophysical contexts, from stellar evolution 
to the structure of our Galaxy. Most recently, interest in 
M dwarfs has been increased further by planet search pro- 
grams. Planets induce higher reflex velocities and deeper 
transits when they orbit and transit M dwarfs rather than 
larger FGK stars, and the habitable zone of the less lumi- 



* Based on observations collected with the FEROS spec- 
trograph at la Silla observatory under ESO programs 
073.D-0802(A), 074.D-0670(A), 078.D-0760(A), and with the 
ELODIE and SOPHIE spectrographs at the Observatoire de 
Haute Provence. 



nous M dwarfs are closer in. Lower mass, smaller, and pos- 
sibly habitable planets are therefore easier to find around 
M dwarfs, and are indeed detected at an increasing pace 
(e.g. |Udry et al.||2007l [Mayor et al.][2009 1 . 

Interesting statistical correlations between the char- 
acteristics of exoplanets and the properties of their host 
stars have emerged from the growing sample of exoplan- 



etary systems (e.g. Endl et al. 2006 Johnson et al. 2007 



Udry fc Santos] |2007[ Bonfils et al. 2011 in prep.). Of 

those, the planet-metallicity correlation was first iden- 
tified and remains the best established: a higher metal 
content increases, on average, the probability that a star 



hosts Jovian planets (Gonzalez 1997 Santos et al. 2001 



2004 Fischer & Valenti 2005). Within the core- accretion 



paradigm for planetary formation, that correlation reflects 
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the higher mass of solid material available to form pro- 
toplanetary cores in the protoplanetary disks of higher 
metallicity stars. The correlation is then expected to ex- 
tend to, and perhaps be reinforced in, the cooler M dwarfs. 
To counterbalance the lower overall mass of their proto- 
planetary disks, those disks need a higher fraction of re- 
fractory material to form similar populations of the proto- 
planetary core. Whether the planet-metallicity correlation 
that seems to vanish for Neptunes and lower mass planets 
around FGK stars ( |Sousa et al.||2008[ [Bouchy et al.|200"9 ) 
persists for Neptune-mass planets around M dwarfs is still 
an open question. 

Our derivation of the first photometric metallicity cal- 
ibration for M dwarfs ( Bonfils et al. 2005 1 was largely 



Valenti et al. 1998 Bean et al. 2006). In this context, 



less direct techniques have been developed to evaluate the 
metal content of M dwarfs. Of those, the most success- 
ful leverage the photometric effects of the very molecular 
bands that complicate spectroscopic analyses. Increased 
TiO and VO abundances in metal-rich M dwarfs shift ra- 
diative flux from the visible range, where these species 
dominate the opacities, to the near infrared. For a fixed 
mass, an increased metallicity also reduces the bolometric 
luminosity. Those two effects of metallicity work together 
in the visible, but, in the [Fe/H] and Te// range of in- 
terest here, they largely cancel out in the near-infrared. 
As a result, the absolute V magnitude on an M dwarf 
is very sensitive to its metallicity, while its near infrared 



motivated by probing their planet-metallicity correlation, 
though only two M-dwarf planetary systems were known 
at the time. A few planet detections later, a Kolmogorov- 
Smirnov test of the metallicity distributions of M dwarfs 
with and without known planets indicated that they only 
had a ~ 11% probability of being drawn from a single par- 



magnitudes are not (Chabrier & Baraffe 2000 Delfosse 



et al. 2000). Position in a color/absolute magnitude di- 



agram that combines visible and near-infrared bands is 
therefore a sensitive metallicity probe, but one that needs 
external calibration. 



We pioneered that approach in Bonfils et al. (2005) 



ent distribution (Bonfils et al. 20071. With an improved where we anchored the relation on a combination spec- 



metallicity calibration and a larger sample of M dwarf troscopic metallicities of early-M dwarfs from Woolf fc 
planets, Schlaufman & Laughlin (2010) lower the proba- Wallerstein (2005) and metallicities, which we measured 



bility that M-dwarf planetary hosts have the same metal- 
licity distribution as the general M dwarf population to 
~ 6%. This result is in line with expectations for the core 
accretion paradigm, but is only significant at the ~2 a 
level. Both finding planets around additional M dwarfs 
and measuring metallicity more precisely will help char- 
acterize this correlation and the possible lack thereof. Here 
we explore the second avenue. 

Measuring accurate stellar parameters from the op- 
tical spectra of M dwarfs unfortunately is not easy. As 
the abundances of diatomic and triatomic molecules (e.g. 
TiO, VO, H2O, CO) in the photospheric layers increases 
with spectral subtype, their forest of weak lines eventu- 
ally erases the spectral continuum and makes a line-by- 
line spectroscopic analysis difficult for all but the earlier 



for the FGK primaries of binary systems containing a 
widely separated M dwarf component. That calibration, 
in terms of the i^s-band absolute magnitude and the 
V — Kg color, results in a modestly significant disagree- 
ment between the mean metallicity of solar-neighborhood 



early/mid-M dwarfs and FGK dwarfs. Johnson & Apps 
(2009) correctly points out that M and (at least) K dwarfs 
have the same age distribution, since both live longer than 
the age of the universe, and that they are therefore ex- 
pected to have identical metallicity distributions. They 
derived an alternative calibration, anchored in FGK-f-M 



binaries that partly overlap the Bonfils et al. (20051 sam- 



ple, which forces the agreement of the mean metallicities 
of local samples of M and FGK dwarfs. Most recently. 



M subtypes. Woolf & Wallerstein (2005 2006) measured 



Schlaufman & Laughlin (2010) have pointed out the im 



atomic abundances from the high-resolution spectra of 67 
K and M dwarfs through a classical line-by-line analysis, 
but had to restrict their work to the earliest M subtypes 
(Te/ / > 3500 K) and to mostly metal-poor stars (median 
[Fe/II]= —0.89 dex). They find that metallicity correlates 
with CaH and TiO band strengths, but do not offer a 
quantitative calibration. 

Although the recent revision of the solar oxygen abun- 
dance ( )Asplund et al.|2009 Caffau et al.|2011 ) has greatly 
improved the agreement between model atmosphere pre- 
diction and spectra of M dwarfs observed at low-to- 
medium resolution (Allard et al. 2010), many visual-to- 



red spectral features still correspond to molecular bands 
that are missing or incompletely described in the opac- 
ity databases that underly the atmospheric models. At 
high spectral resolution, many individual molecular lines 
in synthetic spectra are additionally displaced from their 
actual position. Spectral synthesis, as well, has therefore 
had limited success in analyzing M dwarf spectra (e.g. 



portance of kinematically matching the M and GK sam- 
ples before comparing their metallicity distributions, and 
used stellar structure models of M dwarfs to guide their 
choice of a more effective parametrization of position in 
the Mk^ vs V — Kg diagram. The difference between the 
three calibrations varies slightly across the Herzprung- 
Russell diagram but, on average, the [Johnson fc Apps 
(2009]) calibration is 0.2 dex more metal-rich than Bonfils 
et al j( |2005[ ), and [Schlaufman fc Laughlin] pOlO] is half- 
way between those two extremes. Those discrepancies are 
largely irrelevant when comparing M dwarfs with metal- 
licities consistently measured on any of these three scales, 
but they are uncomfortably large in any comparison with 
external information. 

We set out here to test those three calibrations. For 
that purpose, we have assembled a sample of 23 M 
dwarfs with accurate photometry, parallaxes, and metal- 
licity measured from a hotter companion (Sect. [2]). We 
then perform statistical tests of the three calibrations 
in Sect. |3j and in Sect. [4] we discuss those results and 
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slightly refine the Schlaufman & Laughlin 



(|2010| calibra- 
tion, which we find works best. Section 5| presents our 
conclusions, and an appendix compares our preferred cal- 
ibration against metallicities obtained with independent 
techniques. 

2. Sample and observations 

We adopt the now well-established route of measuring 
the metal content of the primaries of FGK-I-M binaries 
through classical spectroscopic methods, by assuming that 
it applies to the M secondaries. We searched for such bi- 
naries in the third edition of the catalog of nearby stars 



( Gliese & Jahreifi 1991 1, the catalog of nearby wide binary 
and multiple systems ( Poveda et al.||1994 ), the catalog of 
common proper-motion companions to Hipparcos stars 
(Gould & Chaname 2004), and the catalog of disk and 



halo binaries from the revised Luyten catalog (Chaname 
fc Gould|2004l ). To ensure uncontaminated measurements 



of the fainter M secondaries, we required separations of at 
least 5 arcsec. That initial selection identified almost 300 
binaries. We eliminated known fast rotators, spectroscopic 
binaries, pairs without a demonstrated common proper 
motion, as well as systems that do not figure in the re- 
vised Hipparcos catalog ( van Leeuwen|[2"007 ) from which 
we obtained the parallaxes of the primaries. With very few 
exceptions, the secondaries have good JHKs photometry 
in the 2MASS catalog ([Skrutskie et al.||2006|), which we 



therefore adopt as our source of near-infrared photometry. 
The only exception is Gl 551 (Proxima Centauri), which 
has saturated K,, 2MASS measurements and for which 



RIJH photometry was used to calculate metallicity from 



the Casagrande et al. (2008) calibration, as discussed in 
the Appendix. Our final sample contains 23 systems, of 
which 19 have M-dwarf secondaries and four have K7/K8 
secondaries. 

We either measured the metallicity of the primaries 
from high-resolution spectra or adopted measurements 



from the literature which are on the same metallicity scale. 
We obtained spectra for nine stars with the FEROS spec- 



trograph (Kaufer fc Pasquini 1998 ) on the 2.2m ESO/MPI 



telescope at La Silla. We used the ARES program ( Sousa 



et al. 2007) to automatically measure the equivalent 



widths of the Pel and Fe2 weak hnes (< 200 mA) in 



the Fe line list of Sousa et al. (2008). This list is com 



prised of 263 Fe 1 and 36 Fe 2 stable lines, ranging, in 
wavelength, from 4500 to 6890 A. Then, we followed the 
procedure described in Santos et al. (2004): [Fe/H] and 
the stellar parameters are determined by imposing exci- 
tation and ionization equilibrium, using the 2002 version 



of the MOOG (Sneden 1973) spectral synthesis program 



with a grid of ATLAS9 plane-parallel model atmospheres 
( |Kurucz|[T993l ). 



we use the Bessell ( 1991 1 measurements that we trans- 
form into Kg photometry using the equations of [Carpenter] 
(20011). 



Precise optical photometry of the secondaries, to our 
initial surprise, has been less forthcoming, and we sus- 
pected noise in their y-band photometry to contribute 
much of the dispersion seen in previous photometric 
metallicity calibrations. We therefore applied a strict 
threshold in our literature search and only retained pairs 
in which the y-band magnitude of the secondary is mea- 
sured to better than 0.03 magnitude. This criterion turned 
out to severely restrict our sample, and we plan to ob- 
tain y-band photometry for the many systems that meet 
all our other requirements, including the availability of a 
good high-resolution spectrum of the primary. |Mermilliod| 
et al. (1997) has been our main source of Johnson- 
Cousins VRI photometry. For ten sources RI photometry 
was in Weistrop and Kron systems instead of Johnson- 
Cousins. We therefore applied transformations following 
Weistrop] (11975]) and ]Leggett] (]1992|, respectively. The 



For three stars, we used spectra gathered with the 
CORALIE (jQuelozetaLpOOOl spectrograph, on the Swiss 



Euler 1.2 m telescope at la Silla, and SOPHIE (Bouchy & 



The Sophie Team 2006 1 spectrograph, on the Observatoire 
de Haute Provence 1.93 m telescope. For those three stars, 
we use metallicities derived from a calibration of the equiv- 
alent width of the cross correlation function (CCF) of their 
spectra with numerical templates ( Santos et al.|2002[ ) . We 
adopted that approach, rather than a standard spectro- 
scopic analysis, because those observations were obtained 
with a ThAr lamp illuminating the second fiber of the 
spectrographs for highest radial velocity precision. The 
contamination of the stellar spectra by scattered ThAr 
light would affect stellar parameters measured through 
a classical spectral analysis, but is absorbed (to first or- 
der) into the calibration of the CCF equivalent width to a 
metallicity. That calibration is anchored onto abundances 



derived with the Santos et al. (2004) procedures, and has 



been verified to be on the same scale to within 0.01 dex 
( ]Sousa et al.|]20TI] ). 



We adopt 10 [Fe/H] determinations from previous pub- 



2008) 



(Bonfils et al. 


2005 


I the 


Santos et al. 



Sousa et al. 



ods. Finally, we take one metallicity value from |Valenti| 
& Fischer (2005). That reference derived its metallicities 



through full spectral synthesis, and Sousa et al. (20081 
found that they are on the same scale as Santos et al.| 
(J2004]). 



Table [T] lists the adopted stellar parameters (effec- 
tive temperature, surface gravity, micro-turbulence, and 
metallicity) from high-resolution spectra of the primaries. 
Table ]2] lists parallaxes and photometry for the full sam- 
ple, along with their respective references. Columns 1 and 
3 display the names of the primary and secondary stars, 
while columns 2 and 4 display their respective spectral 
types. Column 5 lists the Hipparcos parallaxes of the pri- 
maries with their associated standard errors. Columns 6 
to 11 contain the V{RI)cJHKs photometry of the sec- 
ondary and their associated errors. Column 12 contains 
the bibliographic references for the photometry. 
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Table 1. Stellar parameters measured from the primaries, with the [Fe/H] of the M dwarf secondary inferred from 
the primary. 



Primary 


Secondary 


Teff 

[K] 






log 
[cm s 




] 




6 

[km s 






[Fc/H] 




[Fe/H] 
source 


source 


G153.1A 


G153.1B 


4705 ± 


131 


4 


33 


± 





26 





76 


± 





25 





07 


± 





12 




B05 


G156.3A 


G156.3B 


5394 ± 


47 

























00 


± 





10 


COR 


S08CAL 


G181.1A 


G181.1B 


5332 ± 


22 


3 


90 


± 





03 





99 


± 





02 





08 


± 





02 




SOS 


GllOOA 


GllOOC 


4804 ± 


81 


4 


82 


± 





24 


1 


25 


± 





24 


-0 


28 


± 





03 




New 


G1105A 


G1105B 


4910 ± 


65 


4 


55 


± 





14 





77 


± 





18 


-0 


19 


± 





04 




New 


G1140.1A 


G1140.1B 


4671 ± 


65 


4 


31 


± 





15 





54 


± 





31 


-0 


41 


± 





04 




SOS 


G1157A 


G1157B 


4854 ± 


71 


4 


75 


± 





19 


1 


31 


± 





20 


-0 


16 


± 





03 




New 


G1173.1A 


G1173.1B 


4888 ± 


72 


4 


72 


± 





16 





97 


± 





21 


-0 


34 


± 





03 




New 


G1211 


G1212 


5293 ± 


109 


4 


50 


± 





21 





79 


± 





17 





04 


± 





11 




B05 


G1231.1A 


G1231.1B 


5951 ± 


14 


4 


40 


± 





03 


1 


19 


± 





01 


-0 


01 


± 





01 




New 


G1250A 


G1250B 


4670 ± 


80 


4 


41 


± 





16 





70 


± 





19 


-0 


15 


± 





09 




BOS 


G1297.2A 


G1297.2B 


6461 ± 


14 


4 


65 


± 





02 


1 


74 


± 





01 





03 


± 





05 




New 


G1324A 


G1324B 


5283 ± 


59 


4 


36 


± 





11 





87 


± 





08 





32 


± 





07 




B05 


G1559A 


G1551 


5857 ± 


24 


4 


38 


± 





04 


1 


19 


± 





03 





23 


± 





02 




New 


G1611A 


G1611B 


5214 ± 


44 


4 


71 


± 





06 












-0 


69 


± 





03 




SPO 


G1653 


G1654 


4723 ± 


89 


4 


41 


± 





24 





52 


± 





31 


-0 


62 


± 





04 




SOS 


G1666A 


G1666B 


5274 ± 


26 


4 


47 


± 





04 





74 


± 





05 


-0 


34 


± 





02 




New 


G1783.2A 


G1783.2B 


5094 ± 


66 


4 


31 


± 





13 





30 


± 





19 


-0 


16 


± 





08 




B05 


G1797A 


G1797B 


5889 ± 


32 


4 


59 


± 





06 


1 


01 


± 





06 


-0 


07 


± 





04 




B05 


GJ3091A 


GJ3092B 


4971 ± 


79 


4 


48 


± 





15 





81 


± 





22 





02 


± 





04 




SOS 


GJ3194A 


GJ3195B 


5860 ± 


47 

























00 


± 





10 


SOP 


S08CAL 


GJ3627A 


GJ3628B 


5013 ± 


47 






















-0 


04 


± 





10 


SOP 


S08CAL 


NLTT34353 


NLTT34357 


5489 ± 


19 


4 


46 


± 





03 





91 


± 





03 


-0 


18 


± 





01 




New 



References. [B05] [Bonfils et al.|(|2005[); [C OR] CCF [Fe/H] derived from spectra of the COR ALIE Spectrograph ; [SOSC AL] 
Te// calibration from |Sousa et al.| ( |200S| ); [SOS] |Sousa et al.| ( |200S|); [New] This paper; [SPO] [Valenti fc Fischer] ( [20051 ) ; POP] 
CCF [Fe/H] taken from spectra of the SOPHIE Spectrograph ( [Bouchy fc The Sophie Team|20"06f ! 



3. Evaluating the photometric metallicity 
calibrations 

To assess the three alternative photometric calibrations, 
we evaluated the mean and the dispersion of the difference 
between the spectroscopic metallicities of the primaries 
and the metallicities that each calibration predicts for the 



M dwarf components. As in previous works ( Schlaufman & 



Laughlin 2010 Rojas-Ayala et al. 2010 ), we also computed 



the residual mean square R MSp and the s quared multiple 
correlation coefficient R^p ( Hocking||l976 ). 

The residual mean square RMSp is defined as 



RMSp ^ 

n — p 



SSEp = ^ {y. 



odel 



(1) 



where SSEp is the sum of squared residuals for a p-term 
model, n the number of data points, and p the number of 
free parameters of the model. The squared multiple cor- 



relation coefficient R^^ is defined as 



Rlp^l-{n- 1) 



RMSp 
SST 



ssT = (y^ - 



(2) 



A low RMSp means that the model describes the data 
well, while R^^ close to 1 signifies that the tested model 
explains most of the variance of the data. The R^^ can 
take negative values, when the model under test increases 
the variance over a null model. 



We recall that p should be set to the number of ad- 
justed parameters when a model is adjusted, but instead 
is zero when a preexisting model is evaluated against in- 
dependent data. We are, somewhat uncomfortably, in an 
intermediate situation, with 11, 2, and 12 binary sys- 
tems in common with the samples that define the calibra- 



tions of Bonfils et al. 


(200£ 


>l, 


Johnson & Apps 


(2009 


), and 


Schlaufman & Laughlin 


(■ 


2010 


), and some measurements 



for those systems in common. Our sample therefore is not 
fully independent, and in full rigor p should take some ef- 
fective value between zero and the number of parameters 
in the model. Fortunately, that number, 2 for all three cal- 
ibrations, is a small fraction of the sample size, 23. The 
choice of any effective p between and 2 therefore has lit- 
tle impact on the outcome. We present results for p = 0, 
except when adjusting an update of the Schlaufman & 



Laughlin (2010) calibration to the full sample, where we 
use p = 2 as we should. 

We evaluate the uncertainties on the offset, dispersion, 
RMSp, and R^^ through bootstrap resampling. We gen- 
erated 100,000 virtual samples with the size of our ob- 
served sample by random drawing elements of our sample, 
with repetition. We computed the described parameters 
for each virtual sample, and used their standard deviation 
to estimate the uncertainties. 

Table |3] displays the defining equations of the various 
calibrations, their mean offset for our sample, the dis- 
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persion around the mean value (rms), the residual mean 
square {RMSp), the square of the multiple correlation co- 
efficient (-R^p), as well as their uncertainties. The Mk from 
the B05 calibration is the absolute magnitude calculated 
with the Ks photometric magnitudes and the Hipparcos 
parallaxes. The AM from the B05(2) calibration is the dif- 
ference between the V- and the K-hand mass-luminosity 
relations of Delfosse et al. (2000). In the JA09 calibra- 



0.07 dejQ than 1000 FGK stars, w ith a modest s ignifi 



tion, the AMk is the difference between the mean value 



of [Fe/H] of the main sequence FGK stars from the Valenti 
& Fischer ( 2005[ ) catalog (as defined by a fifth-order poly- 
nomial MS = J2a,{V-Ks)\ where a = {-9.58933, 
17.3952, -8.88365, 2.22598, -0.258854, 0.0113399}), and 
the absolute magnitude in the Ks band. The A(V — Kg) 
in the SLIO and 'This paper' calibrations is the difference 
between the observed V — Kg color and the fifth-order 
polynomial function of Mk^ adapted from the previously 



mentioned formula from Johnson & Apps ( 2009 1 . In this 



case, the coefficients of the polynomial are, in increasing 
order: (51.1413, -39.3756, 12.2862, -1.83916, 0.134266, 
-0.00382023). 

Figure [l] d epicts the different [Fe/H] calibrations from 
Bonfils et aL| (2005[) (a and b), [Johnson fc Ap'^ ((2009) 
(c), Schlaufman & Laughlin (2010) (d), and the calibra- 
tion determined in this paper (e). Table |4] displays the 
metallicity values from spectroscopy and the different cal- 
ibrations, where the individual values for each star can be 
compared directly. 

The bootstrap uncertainties of the parameters (Table 
[3]) show that the rms values are the most robust. The R^p 
parameter, in contrast, has large uncertainties. With our 
small sample size, it therefore does not provide an effective 
diagnostic of the alternative models. 



4. The latest metallicity measurements and 
calibrations 

In this section we discuss the three photometric metal- 
licity calibrations in turn, and examine their agreement 
with our spectroscopic sample. Figure |2] plots the [Fe/H] 
obtained from each calibration against the spectroscopic 
[Fe/H], and it guides us through that discussion. 



4.1. Bonfils et al. (2005) calibration 



As recalled in the introduction, B05 first calibrated po- 
sition in a {(V — i^s) — Mk^} color-magnitude diagram 
into a useful metallicity indicator. That calibration is 
anchored, on the one hand, in spectroscopic metallicity 
measurements of early metal-poor M-dwarfs by [Woolf fc| 
Wallerstein (2005), and on the other hand, in later and 



more metal-rich M dwarfs which belong in multiple sys- 
tems for which B05 measured the metallicity of a hot- 
ter component. The B05 calibration has a ^0.2 dex dis- 
persion. Then, they used the calibration to measure the 
metallicity distribution of a volume-limited sample of 47 
M dwarfs, which they found to be more metal-poor (by 



cance of 2.6 a. As mentioned above, Bonfils et al. (20071 



used that calibration to compare M dwarfs with and with- 
out planets, and found that planet hosts are marginally 
metal-rich. 

For our sample, the B05 calibration is offset by — 0.04± 
0.04 dex and has a dispersion of 0.20±0.02 dex. The nega- 



tive offset is in line with SLIO finding (see Section 4.3 ) that 



B05 generally underestimates the true [Fe/H]. Correcting 
from this —0.04 offset almost eliminates the metallicity 
difference between local M dwarfs and FGK stars. 

SLIO also report that the B05 calibration has a very 
poor under 0.05, and that their own model explains 
almost an order of magnitude more of the variance of their 
calibration sample. In Sect. 3, we noted, however, that R^p 
is a noisy diagnostic for small samples. 

In addition to their more commonly used calibra- 
tion, B05 provide an alternative formulation for [Fe/H]. 
That second expression, labeled B05(2) in Table [3) works 
from the difference between the V- and iC,-band mass- 



luminosity relations of Delfosse et al. (2000). The two B05 



formulations perform essentially equally for our sample, 
with B05(2) having a marginally higher dispersion. In the 
remainder of this paper we therefore no longer discuss 
B05(2). 



4.2. \Johnson & App^ ^200^ calibration 



Johnson & Apps (2009) argue that local M and FGK 



dwarfs should have the same metallicity distribution, and 
accordingly chose to fix their mean M dwarf metallicity to 
the value (—0.05 dex) for a volume-limited sample of FGK 
dwarfs from thelValenti & Fischer ( 2005 ) sample. They de- 



fined a sequence representative of average M dwarfs in the 
{{V — Ks) — Mk^} color-magnitude diagram, and used the 
distance to that main sequence along Mk^ as a metallicity 
diagnostic. They note that the inhomogeneous calibration 
sample of B05 is a potential source of systematics, and 
consequently chose to calibrate their scale from the metal- 
licities of just six metal-rich M dwarfs in multiple systems 
with FGK primary components. 

JA09 present two observational arguments for fixing 
the mean M dwarf metallicity. They first measured [Fe/H] 
for 109 G0-K2 stars (4900<Te//<5900 K) and found 
no significant metallicity gradient over this temperature 
range, from which they conclude that no difference is 
to be expected for the cooler M dwarfs. We note, how- 
ever, that a linear fit to their G0-K2 data set ([Fe/H] = 
9.74 X 10-5 (Teff - 5777) - 0.04) allows for a wide metal- 
licity range when extrapolated to the cooler M dwarfs 
(2700<Teff<3750, for M7 to MO spectral type, with [Fe/H] 
= —0.24 allowed at the 1 a level for an MO dwarf and 
significantly lower than the [Fe/H] offset in B05. More 
importantly, they measured a large (0.32 dex) offset be- 
tween the B05 metallicities of six metal-rich M dwarfs in 



erroneously quoted as a 0.09 dex difference in Johnson & 



Apps (20091 
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Calibration Source + equation 




offset 


rms 


RMSp 








[dex] 


[dex] 


[dex] 




BOS : [Fe/H] = 0.196 - 1.527Mk + 0.091M^ + 1.886(y - K,) - 0.142(V - 




-0.04 ± 0.04 


0.20 ± 0.02 


0.04 ± 0.01 


0.31 ± 0.22 


B05(2) : [Fe/H] = -0.149 - 6.508AM, AM = Massv - MassK 




-0.05 ± 0.04 


0.22 ± 0.02 


0.05 ± 0.01 


0.21 ± 0.34 


JA09 : [Fe/H] = 0.56AMjf - 0.05, AMjf = MS - Mk 




0.14 ± 0.04 


0.24 ± 0.04 


0.06 ± 0.02 


0.03 ± 0.51 


SLIO : [Fe/H] = 0.79A(V - AT,) - 0.17, A(V - K^) = {V - K^)abs - {V - 




0.02 ± 0.04 


0.19 ± 0.03 


0.04 ± 0.01 


0.41 ± 0.29 


This paper : [Fe/H] = 0.57A(y - K^) - 0.17 




0.00 ± 0.04 


0.17 ± 0.03 


0.03 ± 0.01 


0.43 ± 0.23 



Table 4. Spectroscopic metallicity of the primaries and metallicities predicted for the secondary by the different 
calibrations. 



Primary 


Secondary 


Spectroscopic 


[Fe/H] [dex] 
B05 B05(2) 


JA09 


SLIO 


This paper 


G153.1A 


G153.1B 


0.07 


-0.21 


-0.19 


-0.05 


-0.17 


-0.17 


G156.3A 


G156.3B 


0.00 


-0.34 


-0.42 


-0.07 


-0.21 


-0.20 


G181.1A 


G181.1B 


0.08 


-0.22 


-0.30 


0.02 


-0.10 


-0.12 


GllOOA 


GllOOC 


-0.28 


-0.39 


-0.38 


-0.31 


-0.41 


-0.34 


G1105A 


G1105B 


-0.19 


-0.18 


-0.18 


-0.03 


-0.15 


-0.15 


G1140.1A 


GI140.1B 


-0.41 


-0.38 


-0.44 


-0.12 


-0.25 


-0.23 


G1157A 


G1157B 


-0.16 


0.04 


0.13 


0.36 


0.20 


0.10 


G1173.1A 


G1173.1B 


-0.34 


-0.27 


-0.25 


-0.14 


-0.25 


-0.23 


G1211 


G1212 


0.04 


-0.08 


-0.09 


0.15 


0.04 


-0.02 


G1231.1A 


GI231.1B 


-0.01 


-0.11 


-0.06 


0.15 


0.01 


-0.04 


G1250A 


G1250B 


-0.15 


-0.18 


-0.14 


0.04 


-0.09 


-0.11 


G1297.2A 


GI297.2B 


0.03 


0.00 


0.05 


0.27 


0.13 


0.05 


G1324A 


G1324B 


0.32 


-0.01 


0.04 


0.34 


0.22 


0.11 


G1559A 


G1551 


0.23 


0.06 


-0.08 


0.19 


0.20 


0.10 


G1611A 


G1611B 


-0.69 


-0.30 


-0.40 


-0.64 


-0.81 


-0.64 


G1653 


G1654 


-0.62 


-0.27 


-0.26 


-0.07 


-0.19 


-0.18 


G1666A 


G1666B 


-0.34 


-0.09 


-0.14 


0.12 


0.02 


-0.03 


G1783.2A 


GI783.2B 


-0.16 


-0.15 


-0.15 


0.02 


-0.10 


-0.12 


G1797A 


G1797B 


-0.07 


-0.02 


0.04 


0.25 


0.11 


0.03 


GJ3091A 


GJ3092B 


0.02 


-0.15 


-0.22 


-0.15 


-0.27 


-0.25 


GJ3194A 


GJ3195B 


0.00 


-0.19 


-0.14 


0.04 


-0.10 


-0.12 


GJ3627A 


GJ3628B 


-0.04 


-0.10 


-0.06 


0.16 


0.03 


-0.03 


NLTT34353 


NLTT34357 


-0.18 


-0.34 


-0.38 


-0.10 


-0.22 


-0.21 



multiple systems and the spectroscopic metallicities which 
they measured for their primaries. This robustly points 
to a systematic offset in the B05 calibration for metal- 
rich M dwarfs, but does not directly probe the rest of 
the (Te//, [Fe/H]) space. We do find that the JA09 cali- 
bration is a good metallicity predictor for our sample at 
high metallicities, where its calibrator was chosen. With 
decreasing metallicity, that calibration increasingly over- 
estimates the metallicity, however, as previously pointed 
out by SLIO (see below). Quantitatively, we measure a 
-|-0.14± 0.04 dex offset for our sample and a dispersion of 
0.24 ±0.04. 



4.3. Schlaufman & Laugh lin (2010^ calibration 



Schlaufman & Laughlin (20101 improve upon B05 and 



JA09 in two ways. They first point out that, for M and 
FGK dwarfs to share the same mean metallicity, matched 
kinematics is as important as volume completeness. Since 



the various kinematic populations of our Galaxy have very 
different mean metallicities, the mean metallicity of small 
samples fluctuates very significantly with their small num- 
ber of stars from the metal-poor populations. To over- 
come this statistical noise, they draw from the Geneva- 
Copenhagen Survey volume-limited sample of F and G 
stars a subsample that kinematically matches the volume 
limited sample of M dwarfs used by JA09. They find a 
~ —0.14 ± 0.06 dex mean metallicity for that sample, 
0.09 dex lower than adopted by JA09. However, they only 
used that sample to verify that the mean metallicity of M 
dwarfs in the solar neighborhood is well defined. In the 
end, the M dwarfs within a sample of binaries with an 
FGK primary that they used to fix their calibration are 
not volume-limited or kinematically-matched, but their 
mean metallicity ([Fe/H] — —0.17 ± 0.07) is statistically 
indistinguishable from the mean metallicity of the volume- 
limited and kinematically-matched sample. 
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B05 calibration 
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JA09 calibration 
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Fig. 1. The different [Fe/H] calibrations from 



Bonfils et al. 



(20051 (a and b), Johnson fc Apps (2009) (c), 



Schlaufman 



& Laughlin (2010) (d), and the calibration determined in this paper (e). In each upper panel, the blue/black dots 



represent the data points. The black line depicts a fit to the data except in panel (a), where the calibrated [Fe/H] is 
shown as isometallicity contours. The lower subpanels show the difference between the calibrated and the spectroscopic 
metallicity. The black dashed lines represent the null value, and the red dotted line represents the mean difference for 
that calibration. 



Second, they use stellar evolution models to guide 
their parametrization of the color-magnitude space. Using 
[Fe/H] isocontours for the [Baraffe et alT (1998) models, 
they show that in a {(^ — Kg) — Mk^ } diagram, changing 
[Fe/H] affects {V — Kg) at an essentially constant Mk,- 
The metallicity is therefore best parametrized by (V—Ks), 
and their calibration uses a linear function of the {V — Kg) 



distance from a nominal sequence in the {{V—Ks) — Mk^} 
diagram. They do not force any specific mean metallicity, 
but verify a posteriori that it matches expectations. 

We measure a 0.14 ± 0.02 dex dispersion for the SLIO 
sample against their calibration, but that calibration has a 
significantly higher dispersion of 0.19±0.03 for our valida- 
tion sample. That increased dispersion reflects our sample 
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Fig. 2. [Fe/H] estimated from the the cahbrations versus spectroscopic metalhcity.The blue dots with error bars 
represent the data points. The black line depicts a one-to-one relationship. 



probing a wider metallicity range than SLIO, as verified 
by computing the dispersion of an 18 star subsample that 
matches the metallicity range of the SLIO sample. That 
dispersion is 0.15 ± 0.03 dex, and indistinguishable from 
0.14 ± 0.02 dex for the SLID sample. The increased dis- 
persion for a wider metallicity range suggests that a linear 
function of {V — Kg) does not fully describe metallicity. 
We also measure an offset of 0.02 ± 0.04 dex. Offset and 
rms both improve over either of the BOS and JAlO cali- 
brations. 



expression is adopted from Schlaufman & Laughlin (2010 ), 
who adapted an Mk, vs {V — Ks) formula from Johnson 



4.4. Refining the \Schlaufman & Laugh lin\ (20 IC^ 
calibration 

We produced updated coefficients for the SLIO prescrip- 
tion, using the RMSp free parameter p = 2 (see Sect. [s]). 
The expression for the new metallicity calibration is 



fc Apps| ( |2009[ ). 

Table |3| shows limited differences between this new fit 
and the original SLIO calibration. The dispersion of the 
new fit is tighter by just 0.02 dex (0.17±0.03 dex instead of 
0.19 ±0.03), and the offset is now 0.00 ±0.04, as expected. 
The Rip value is similar (0.43 ± 0.23 vs 0.41 ± 0.29) and 
uncertain. Readjusting the coefhcients therefore produces 
a marginal improvement at best. 

The dispersion, in all panels of Fig. [T] is well above the 
measurement uncertainties. Those therefore contribute 
negligibly to the overall dispersion, which must be domi- 
nated by other sources. 

As can be seen in Fig. [2) B05 or B05(2) tend to under- 
estimate [Fe/H], while the JA09 calibration clearly over- 
estimates [Fe/H] except at the highest metallicities. 



[Fe/H] = 0.57A(y - K^) - 0.17, 



(3) 



A{V - Ks) ^{V-Ks 



[V-Ks 



where {V — Ks)obs is the observed V — Ks color and 
{y — Ks)iso is a fifth-order polynomial function of Mfc^ 
that describes the mean main sequence of the solar neigh- 



borhood from the Valenti & Fischer ( 2005 1 catalog. This 



5. Summary 

We have assembled a sample of M dwarf companions to 
hotter FGK stars, where the system has an accurate par- 
allax and the M dwarf component has accurate V and Kg- 
band photometry. Using the metallicities of the primaries, 
newly measured or retrieved from the literature, and the 
assumption that the two components have identical initial 
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compositions, we compared the dispersion s of the|Bonfils 
et al.|(|2005|),|Johnson & Appsl (|2009|), and ISchlaufman & 



Laughlin (2010) photometric metallicity calibrations. We 



find that the Schlaufman & Laughlin ( 2010 ) scale, which is 



intermediate between Bonfils et al. ( 2005 ) and Johnson & 



Apps (2009), has the lowest dispersion. We slightly refine 



that relation, by readjusting its coefficients on our sample. 

We find that our tight selection of binaries with accu- 
rate parallaxes and photometry sample has insignificantly 
reduced the dispersion of the measurements around the 
calibration compared to looser criteria. This suggests that 
the dispersion, hence the random errors of the calibration, 
is not defined by measurement uncertainties but instead 
reflects intrinsic astrophysical dispersion. Nonlinearities in 
the metallicity dependence of the V — Kg color are likely 
to contribute, as suggested both by atmospheric models 
(Allard, private communication) and by the increased dis- 
persion that we measure over a wider metallicity range. 
They are, however, unlikely to be the sole explanation, 
since we see dispersion even in narrow areas of the color- 
magnitude diagram. Stellar evolution cannot significantly 
contribute, since early-M dwarfs evolve rapidly to the 
main sequence and they remain there for much longer than 
a Hubble time, but rotation and magnetic activity could 
play a role. Unless, or until, we develop a quantitative 
understanding of this astrophysical dispersion, the photo- 
metric calibration approach may therefore have reached an 
intrinsic limit. Those calibrations also have the very prac- 
tical inconvenience of needing an accurate parallax. This 
limits their use to the close solar neighborhood, at least 
until the GAIA catalog becomes available in a decade. 

Alternative probes of the metallicities of M dwarfs are 
therefore obviously desirable. One obvious avenue is to 
work from higher spectral resolution information and to 
identify spectral elements that are most sensitive to metal- 
licity and others that are most sensitive to effective tem- 
perature. We are pursuing this approach at visible wave- 



lengths (Neves et al. in prep.), as do Rojas-Ayala et al. 
( 2010[ see Appendix A. 2 1 in the near infrared, with en- 
couraging results in both cases. 
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Appendix A: Other methods 



A.l. Calibration of Casagrande et al. (2008) 



In Sect. |4]we described the photometric metallicity cali- 
brations in detail. Casagrande et al. (2008) devised a com- 
pletely different technique, based on their previous study 



of FGK stars using the infrared flux method ( Casagrande 



et al. 2006), to determine effective temperatures and 



metallicities. The infrared flux method uses multiple pho- 
tometry bands to derive effective temperatures, bolomet- 
ric luminosities, and angular diameters. The basic idea of 



IRFM (Blackwell & Shallis 1977) is to compare the ratio 



between the bolometric flux and the infrared monochro- 
matic flux, both measured on Earth, to the ratio be- 
tween the surface bolometric flux (cx aTeff^) and the 
surface infrared monochromatic flux for a model of the 
star. To adapt this method to M dwarfs, |Casagrande| 
et al. (2008) added optical bands, creating the so-called 



MOITE, Multiple Optical and Infrared TEchnique. This 
method provides sensitive indicators of both temperature 
and metallicity. The proposed effective temperature scale 
extends down to 2100-2200 K, into the L-dwarf limit, and 
is supported by interferometric angular diameters above 
^ 3000K. Casagrande et al. ( |2008 ) obtain metallicities by 
computing the effective temperature of the star for each 
color band {V{RI)cJHKs) for several trial metallicities, 
between —2.1 and 0.4 in 0.1 dex steps, and by selecting the 
metallicity that minimizes the scatter among the six trial 



effective temperatures. Casagrande et al. ( 2008 ) estimate 



that their total metallicity uncertainty is 0.2 to 0.3 dex. 

The MOITE method does not reduce into a closed 
form that can be readily applied by others, but Luca 
Casagrande kindly computed MOITE [Fe/H] values for 
our sample (Table |A.l ) . We evaluated the calibration in 
the same manner as in Sect. |3] and obtained a value of 
-0.11±0.07 dex for the offset, 0.32±0.06 dex for the rms, 
0.10±0.04 dex for the RMSp, and -1.09±1.45 for the R^p. 
From these values and from Fig. |A.1[ we can observe that 
the Casagrande et al. ([2008 ) calibration has a higher rms 
and RMSp and a poorer R^^ than the three photometric 
calibrations, consistently with the high metallicity uncer- 



tainty referred by Casagrande et al. ( 2008 1 . The negative 



R^p value formally means that this model increases the 
variance over a constant metallicity model, but as usual 
R^ is a noisy diagnostic. 



A. 2. \Rojas-Ayala et al.\ {2010^ calibration 



Rojas-Ayala et al. (2010) have recently published a novel 



and potentially very precise technique for measuring M 
dwarf metallicities. Their technique is based on spectral 
indices measured from moderate-dispersion (i? ~ 2700) 
K-hand spectra, and it needs neither a V magnitude nor 
a parallax, allowing measurement of fainter (or/and far- 
ther) stars. They analyzed 17 M dwarf secondaries with 
an FGK primary, which also served as metallicity cali- 
brators, and measured the equivalent widths of the Nal 



12 



Neves et al.: Metallicity of M dwarfs. II. 



Table A.l. Metallicity values from spectroscopy and ob- 



tained using the method of Casagrande et al. (2008) (COS 
in this table). 



Primary 


Secondary 


[Fe/H] [dex 
Spectroscopic 


COS 


G153.1A 


G153.1B 


0.07 


-0.07 


G156.3A 


G156.3B 


0.00 


-0.21 


G181.1A 


G181.1B 


0.08 


-0.08 


GllOOA 


GllOOC 


-0.28 


-0.10 


G1105A 


G1105B 


-0.19 


-0.30 


G1140.1A 


G1140.1B 


-0.41 


-0.30 


G1157A 


G1157B 


-0.16 


-0.10 


G1173.1A 


G1173.1B 


-0.34 


-0.20 


G1211 


G1212 


0.04 


-0.21 


G1231.1A 


G1231.1B 


-0.01 


-0.28 


G1250A 


G1250B 


-0.15 




G1297.2A 


G1297.2B 


0.03 


0.00 


G1324A 


G1324B 


0.32 


-0.20 


G1559A 


G1551 


0.23 




G1611A 


G1611B 


-0.69 


-0.40 


G1653 


G1654 


-0.62 


-0.30 


G1666A 


G1666B 


-0.34 




G1783.2A 


G1783.2B 


-0.16 


-0.30 


G1797A 


G1797B 


-0.07 


-0.90 


GJ3091A 


GJ3092B 


0.02 


-0.30 


GJ3194A 


GJ3195B 


0.00 


-0.60 


GJ3627A 


GJ3628B 


-0.04 


-0.20 


NLTT34353 


NLTT34357 


-0.18 


0.19 



doublet (2.206 and 2.209 fim), and the Cal triplet (2.261, 
2.263 and 2.265 fim). With these measurements and a 
water absorption spectral index sensitive to stellar tem- 
peratures, they constructed a metallicity scale with an ad- 
justed multiple correlation coefficient greater than the one 
ofl^ 



Schlaufman & Laughlin (2010) {Rip = 0.63), and also 
with a tighter RMSp of 0.02 when compared to other stud- 



ies (0.05, 0.04, and 0.02 for Bonfils et al 



2005 



Johnson 



|fc Apps||2009| and [Schlaufman fc Laughlin||2010| respec- 
tively). The metallicity calibration is valid over -0.5 to 
-1-0.5 dex, with an estimated uncertainty of ±0.15 dex. 



A test of the Rojas-Ayala et al. (2010) calibration 
for our full sample would be very interesting, but is 
not currently possible for lack of near-infrared spec- 
tra for most of the stars. Seven of our stars, how- 
ever, have their metallicities measured in |Rojas-Ayala| 



et al.| ( |2010[ ) (Gl 212, Gl 231. IB, Gl 250B, Gl 324B, 
G1611B, G1783.2B, and Gl 797B with predicted [Fe/H] of 
0.09,-0.05,-0.04,0.30,-0.49,-0.19, and -0.06 dex, re- 
spectively). We find a dispersion of only 0.08 dex and an 
offset of 0.04 dex offset between our spectroscopic mea- 
surements of the primaries and the |Rojas-Ayala et al.| 
( 2010 ) metallicities of the secondaries. These numbers are 



extremely encouraging, but still have little statistical sig- 
nificance. They will need to be bolstered by testing against 
a larger sample and over a wider range of both metallicity 
and effective temperature. 




-0.6 -0.4 -0.2 0.0 0.2 
[Fe/H] SPEC 




-0.6 -0.4 -0.2 0.0 0.2 
[Fe/H] SPEC 



Fig. A.l. [Fe/H] obtained with the Casagrande et al. 



(2008) method versus the spectroscopic metallicity. The 
blue dots with error bars represent the data points. The 
black line depicts a one-to-one relationship. The metal- 
licity difference between the values of the calibrations 
and the spectroscopic measurements is shown below each 
[Fe/H]-[Fe/H] plot. The black dashed line is the zero point 
of the difference, and the red dotted line represents the av- 
erage of the metallicity difference. 



